Kinetic electrostatic structures in current-carrying pair plasmas 
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The existence and properties of phase-space structures in current-carrying pair plasmas is 
studied by means of the finite amplitude expressions of the pseudo-potential method. Em- 
phasis is given to double layers, solitary structures and periodic waves. The results can be 
applied to electron-positron plasmas as well as to plasmas containing heavier charged parti- 
cles with the same mass and opposite charges, such as fullerene ions and dust particles. They 
can also help to understand numerical simulations that showed the spontaneous formation 
of phase-space holes in linearly stable regimes. 
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I. INTRODUCTION 



The investigation of collective phenomena in pair plasmas, consisting of two species with the 
same mass is gaining interest among the plasma physics community. There are several reasons for 
that: the first one is that electron-positron pair plasmas are present in many astrophysical contexts, 
such as the early universe, neutron stars and active galactic nuclei. On the other hand, recent 
experiments with fullerene pair plasmas [1] have allowed the investigators to observe collective 
phenomena in the laboratory, something that was earlier impeded by the short annihilation times 
of electron-positron laboratory plasmas But another explanation for the outstanding 

interest in pair plasmas is their relatively simple theoretical description, which makes them a 
particularly good subject of investigations for questions in fundamental plasma theory. 

This property of pair plasmas has been emphasized in earlier works that investigated the non- 
linear instability and saturation of current-carrying plasmas [^Q]- In these works, use was made of 
the equality of time-scales in the dynamics of both species to accelerate the onset and evolution of 
instability in a linearly stable regime. It was also found that, after a turbulent transient stage, the 
system reached a stable state in which phase-space structures were present coexisting with distri- 
bution functions one of which being characterized by a flat trapped particle distribution. Kinetic 



structures in pair plasmas were on the other hand investigated in Refs. [7( and |8(, two works on 
which the present article relies. However, in the first, analytical one the investigations were limited 
to small amplitude waves and a complete symmetry in the parameters describing the distribution 
functions of both species whereas in the second, numerical one only solitary phase-space holes 
were considered. Both works did not consider current-carrying plasmas. It is thus appropiate to 
remove these restrictions in order to arrive at a more general theory of electrostatic trapping in 
pair plasmas, which is the purpose of the present article. 

This paper is organized as follows. In section [H] we provide a theoretical background about the 
governing equations and the procedure to obtain finite amplitude equilibrium solutions. Sections 
IIIIl IIVI and are dedicated, respectively, to the study of double layers, phase-space holes and 
periodic solutions. Finally, the results and conclusions are summarized in section IVII 

II. THEORETICAL BACKGROUND 

We consider a drifting collisionless pair plasma with equal temperatures for both particle species. 
The dynamics is governed by the ID Vlasov-Poisson system 

[d t + vd x ±d x $(x,t)d v }f T (x,v,t) =0, (la) 



dl x ${x, t) = J dvf-- J dvf + = n^- n + , (lb) 

where space x, time t, velocity v, the distribution functions f± and electric potential have been 
normalized by the Debye length Xd, the inverse plasma frequency w" 1 , the thermal speed Vr, 
uq/Vt and T/e, respectively. Here, Ad = {n^e 2 / 6qT) 1 / 2 , uj p = {n^e 2 / earn) 1 / 2 and Vr = (T/m) 1 / 2 , 
where no is the unperturbed particle density of each species, e is the magnitude of the electron 
charge, m is the mass and T is the temperature (in Joules) of the two species in the unperturbed 
state. The distribution functions of the homogeneous, unperturbed state are in the center-of-mass 
frame / 0t = (l/v^exp {-{v T v D /2) 2 /2). 

Electrostatic structures in plasmas have been throughly investigated by means of the pseudo- 
pqtential method, first introduced in the kinetic regime in Ref. $] and further developed in Befs.liol. 
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12l (see also Refs. QandQ for a review). In this method, the distribution functions for equilibria 



in the wave frame depend on the constants of motion in a prescribed form, namely 
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where 9{z) represents the Heaviside step function, N± are normalization constants, = vd/2 z F v o 
and a and (3 are the trapping parameters of positive and negative ions, respectively. The single 

2 

particle energies, which are constants of motion, are given by := \ — where we defined 
<£_ := <3? and := ^ — <I> respectively, where \& is the maximum value of $. The separatrix in 
the phase space of both species is then given by e T = 0, separating free (e^ > 0) from trapped 
(e T < 0) particles. 

The distribution functions given by © can be integrated in velocity, and yield the particle 
densities as functions of the electrostatic potential as n+(<]?) = N + nQ(v^,a,^> — <£), n_(<£) = 
N-no(v-, f3, $), where 



n (u, (3, $) := exp(-n 2 /2) [F {u 2 /2, $) + T (J3, $)] 



(3) 



The definitions of the special functions F(v 2 /2, $) and T(/3, $), which represent the contribution 
of free and trapped particles, respectively, to the density, are presented in Refs. Note that 

no(u, f3, 0) = 1 for any u, (3. The Poisson equation is now solved by defining the classical potential 
V(<&) that satisfies $"(x) = n_($) - n + ($) =: — 9V($)/9$. Multiplying both sides by <&'(x) and 
integrating once, we have <£'(x) 2 /2 + V(3>) = 0, where the classical potential is 



V($) = N + [V (v+,a, ¥) - V (v+, a, * - $)] - N-V (y-,P, $) 
and we have defined 

Vb(«, /3, *) := exp (-u 2 /2) [P(/3, $) - 1 + (u 2 /2, 0, $)] . 



(4) 



(5) 



We note that V(0) = 0. The expressions for the special functions P{(3, <J>) and H(u 2 /2, 0, <I>) are 
also found in Ref. jj], To find acceptable, physical solutions two conditions have to be imposed 
upon V(&): a) V(<&) <0if0< < I ) < , I / and b) V(^) = 0. The second condition is usually referred 
to as the nonlinear dispersion relation (NDR), as it links the amplitude of the structure to its phase 
speed. 
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FIG. 1: Existence curves of double layers (DL) in a non-drifting pair plasma. The solid line represents 
standing DLs while the dashed one stands for propagating DLs with a phase speed Uq = 1- Note that for 
both cases two branches exist, called symmetric and asymmetric branches. Note also that there are no 
asymmetric solutions in the propagating case if a > 0.55. 



III. DOUBLE LAYERS 



Now we make use of the expressions presented in the previous section to obtain equilibrium 
solutions of the Vlasov-Poisson system. First, we look for double layer (DL) solutions. Double 

Sers are configurations of phase-space which are associated with a monotonic step-like potentials 
,Q. 

In order to have a DL, the densities of both species must be equal as x — > ±oo, which means 
<]? — ► and —* ^. Let us assume that the densities are unity at $ = 0. This implies that iV_ = 1 
and N+ = l/no(v+,a,^). Under this assumption, the condition n_(^) = n + (^) reduces to 

no(«+,«,*)no(«-,/9,*) = 1- (6) 

This equation has to be solved simultaneously with the NDR for the unknowns uq and in order 
to obtain valid solutions. 



A. Non-drifting plasma 

n 

In Ref. |2||, the existence of DLs in pair plasmas was discussed. They were limited, however, to 
the small amplitude limit and with the further assumption of complete symmetry between positive 
and negative particles (a = (5). Here we show that, even if we drop such limitations, we can still 



5 


-5 

-10 
100 



& 1 



0.01 



quasisymmetric 




asymmet r ic 



-1.5 



-1 -0.5 

a 



0.5 



quasisymmetric ^ ..-^J 






:==s vr. - 

asymmetric 


1.1,1, 



FIG. 2: Existence curves of double layers (DL) in a current-carrying pair plasma with = 2. The solid line 
represents standing DLs while the dashed one stands for propagating DL with a phase speed vq = —0.25. 
Note that for both cases two branches exist, which are the degenerated form of the symmetric and asymmetric 
branches in the case vd = 0. Note also that the curves for — 1, vq = are equivalent to that of vd = 0, 
vq = 1, plotted in Fig. (see text). 

find DL solutions. Their existence curves are plotted in Fig. ^ Note that there are two different 
branches, labelled "symmetric" and "asymmetric". 

The symmetric branch represents solutions for which a = (3. It exists only for a < and admits 
DLs of arbitrary strength. In the limit f < 1 (|a| = \/3\ ^> 1) it corresponds to the solutions of 
Ref. 7], which, for vq = satisfy (1 — a) = 3t/tt / 4\/ r $> . For finite amplitudes we observe that the 
relationship between the defining parameters is very well approximated by an expression of the 
form |a| 1 I /7 = C, where for standing DLs (vq = 0) we found 7 ~ 0.59, C ~ 1.08. This implies that 
the DL becomes the stronger the more flat both trapped particle distributions are. 

The asymmetric branch, on the other side, represents solutions for which the complete symmetry 
between the species is broken. This branch contains also solutions for which a > 0. In the special 
case a = 0, the distribution function of positive particles is flat in the trapped range while trapped 
negative particles form a dip. 

The most relevant physical difference between both branches is that, while solutions in the 
symmetric branch do not present a jump in the densities, this is not true for asymmetric solutions. 



B. Current-carrying plasma 



The picture is changed when we look for double layers in a current-carrying plasma, a case which 
was not discussed in Ref. 7]. The parameter curves representing possible solutions for vd = 2 are 
shown in Fig. 

First of all we note that the expressions for rio(u,{3, 4>) and Vo(u,(3, do only depend on 
u 2 . Therefore we can make a complete correspondence between standing structures in a current- 
carrying plasma (v- = u+ = v_d/2) and propagating structures in a non-current-carrying plasma 
(— V- = i>+ = vq) just by interchanging v_d/2 and i«o- This explains that the solid lines in Fig. |2 
are the same as the dashed ones in Fig. Q Note however that, although the location of solutions 
in parameter space are exactly the same, both represent very different kinds of solutions. 

Nevertheless, if we look for propagating structures in current-carrying plasmas, the symmetry is 
broken and the correspondence is no longer valid. An example of the location in parameter space 
of propagating DLs is given by the dashed line in Fig. |2j Note that we still have two different 
branches but now, as vq ^ imposes a further asymmetry, we do not have a family of solutions 
with a = (3. We can still call the two branches of solutions "degenerated symmetric branch" and 
"degenerated asymmetric branch". The degenerated symmetric branch exists only for a < and 
requires increasing amplitudes as a — > 0. Here we find also a power law of the form |a| 7x I/ = C 
where for vd = 2, vq = —0.25 we obtain 7 « 1.015, C m 6.11. 



IV. SOLITARY PHASE-SPACE HOLES 



The existence and properties of solitary holes in non-drifting pair plasmas was already studied 



in Ref. Therefore here we focus on the existence of solitary holes in current-carrying pair 
plasmas. As the symmetry arguments mentioned in the past section hold for holes as well as for 
DLs, standing holes in a current-carrying plasma correspond one-to-one to propagating holes in a 
plasma with vo = 0. 

Without loss of generality, we will consider here holes in /_, corresponding to bell-like electro- 
static potentials, i.e. V'(0) = 0, V'{^) > 0. We will also consider only the case of a flat trapped 
range in /+ (a = 0). There are several reasons to give special attention to these cases: a) there are 
large numerical evidences that such configurations are more stable and therefore potential attrac- 
tors of the dynamics. For example, in Ref. the evolution of two colliding holes was simulated 
which resulted in a partially flat distribution function, and b) recent simulations of nonlinear insta- 
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FIG. 3: The amplitude of the electrostatic potential for holes in /_ as a function of the trapping parameter 
(3 for different phase velocities. The drift between both species was vu = 2 and a flat trapped range is 
assumed in /+ (a = 0). 

bility 5] and the subsequent turbulence in current-carrying pair plasmas show that this turbulence 
decays towards a stable hole equilibrium in which one of the species presents a flat distribution in 
the resonant range. 

In Fig. |21 we present the curves in the /3, parameter space where solitary hole solutions exist 
in a plasma with vd = 2 for different values of the phase velocity vq. No solutions are possible 
with P > 0. Note that even for large amplitudes, the curves can be approximated by a power law 
of the form |/3| 7 \I' = C, which is actually exact for small amplitudes. 




V. PERIODIC STRUCTURES 



Periodic structures can also appear in pair plasmas. This is important because they are excited 
in experiments such as Ref. ij is usually periodical and also because most numerical codes impose 
periodic boundary conditions in space, giving preference to periodic structures. 

In Ref. [7j a dispersion relation was found for harmonic waves that smoothly joined the limits 
of slow acoustic modes uj/y/2 ~ 0.924A: and plasma waves ujj\pl ~ 1. To understand how the 
finite amplitude of the potential affects these results, we will use a = f3 = 1 — Vq, which for small 
amplitudes gives raise to harmonic waves (B = in Ref. Q). 

Figure0]shows the dispersion relation of these waves with different values of the amplitude of the 
electrostatic potential. There we notice that for large amplitude waves the minimum phase velocity, 
corresponding to the slow acoustic mode, is considerably smaller than the one corresponding to 
very small amplitudes (vq = 0.924\/2 when ^ <C 1). This slowing effect is shown in Fig. where 
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FIG. 5: Velocity of the slow acoustic branch in the limit of large wavelengths as a function of the amplitude 
of the electrostatic potential. 



the velocity of the slow acoustic mode is plotted as a function of the amplitude of the electrostatic 
potential. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we investigated the existence of electrostatic structures in pair plasmas with a drift 
between the species due e.g. to an external electric field. This results extend and complement those 



of Refs.[3] and Q] in several aspects: we do not restrict ourselves to the case of perfect symmetry 
between the species (a = /?), we do not base the results in the small amplitude limit <C 1) and, 
finally we do not consider only plasmas without a drift velocity between the species (vo = 0). By 
removing those limitations we have explored a wide range of possible structures. We have focused 
on some remarkable results that may help to interprete future experimental and numerical data. 

As relevant outcome, we emphasize the existence of asymmetric double layers, that exhibit a 
jump in the densities, as well as the dependence between the velocity of the slow acoustic mode 
and the potential amplitude for periodic waves. This latter result can easily be compared with 
experimental data extracted from a laboratory setup like that of Ref. Q], something that we propose 
here and that would yield interesting results. 

We would like also to underline that the results presented here are also connected with the 
studies about the nonlinear instability and saturation of pair plasmas, presented elsewhere [j]]. 
The connection is twofold: a) in that reference it was shown that nonlinear stability is triggered 
by small amplitude phase-space holes and b) the final stable equilibrium is indeed a structured 
equilibrium that should be studied with the tools presented here. 
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